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We show that high energy hadronic reactions which contain a rapidity gap and a 
hard subprocess have a specific dependence on the kinematic variables, which results in 
^ ' a characteristic behaviour of the survival probability of the gap. We incorporate this 

' mechanism in a two-channel eikonal model to make an essentially parameter-free estimate 

of diffractive dijet production at the Tevatron, given the diffractive structure functions 
measured at HERA. The estimates are in surprising agreement with the measurements 
of the CDF collaboration. We briefly discuss the application of the model to other hard 
processes with rapidity gaps. 



1 Introduction 



There has been much interest in the probabihty of rapidity gaps in high energy interactions to 
survive, since they may be populated by secondary particles generated by rescattering processes, 
see, for example, The effect can be described in terms of screening or absorptive 

corrections. To the best of our knowledge, the term survival probability was introduced by 
Bjorken who estimated the probability using 

^ / |.M(s,&)pe-"W £b 

J \M{s,b)\^cm ' ^ ^ 

where A4 is the amplitude (in impact parameter b space) of the particular process of interest 
at centre-of-mass energy ^/s. Q is the opacity (or optical density) of the interaction of the 
incoming hadron^. 

It is perhaps more accurate to use the term "suppression factor" of a hard process ac- 
companied by a rapidity gap, rather than "survival probability". It depends not only on the 
probability of the initial state to survive, but is sensitive to the spatial distribution of par- 
tons inside the incoming hadrons, and thus on the dynamics of the whole diffractive part of 
the scattering matrix. It is important to note that the suppression factor S*^ is not universal, 
but depends on the particular hard subprocess, as well as the kinematical configurations. In 
particular, S'^ depends on the nature of the colour-singlet (Pomeron or W/Z boson or photon) 
exchange which generates the gap as well as on the distributions of partons inside the proton 
in impact parameter space 0, |ll], [1^. In this paper we emphasize the importance of the 



dependence on the characteristic momentum fractions carried by the active partons in the col- 
liding hadrons. This leads to a much richer structure of the probability of rapidity gaps in 
processes mediated by colour-singlet t-channel exchange. The framework was introduced long 
ago0 ||15|, , but only with the advent of rapidity gap events being observed in hard processes 
at the Tevatron and at HERA, is this rich physics now revealing itself. 

In Section 2 we briefly review the general framework, and, in particular, discuss a two- 
channel partonic model of diffraction. Measurements of diffractive dijet production with a 
leading antiproton have been made recently by the CDF collaboration at the Tevatron. 
This is an ideal process with which to compare the specific predictions of the models for high 
energy diffraction. In Section 3 we specify the partonic structure of the diffractive two-channel 
eigenstates. To set the scene for our main study we first, in Section 4, discuss diffractive dijet 
production assuming, for the moment, that rescattering corrections may be neglected. As was 



emphasised in Ref. ^7j, the calculation of the cross section, based on factorization in terms of 
diffractive structure functions obtained from HERA data, indicates a large discrepancy with the 
CDF measurements — both in the normalisation and in the shape of the observed distribution. 



^That is i[l — exp(— J7/2)] is the usual elastic scattering amplitude in impact parameter space, il/2 is 
frequently called the eikonal. 

^Reviews can be found, for example, in Refs. [l^ . 
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The calculation lies about a factor of 10 above the data; the precise discrepancy depends on 
the kinematic domain. In Section 5 we include rescattering corrections. Clearly these will 
decrease the predictions, since now the rapidity gaps may be populated by secondary particles. 
To allow for rescattering we use the two-channel eikonal model, reviewed in Sections 2 and 
3, with parameters previously determined in a global description of the total, elastic and soft 
diffractive data available in the ISR to Tevatron energy range 0. In this way we are able to 
make an essentially parameter-free prediction of both the normalisation and the shape of the 
CDF diffractive dijet data. In Section 6 we discuss the application of the model to other hard 
processes with rapidity gaps, but on a less quantitative level than for dijet production. In all 
cases the specific rescattering corrections are in the direction to improve the description of the 
data. Finally, in Section 7, we present our conclusions. 



2 Inelastic diffraction and diffractive eigenstates 

In order to deduce the behaviour of inelastic diffraction, we start with the s-channel unitarity 
relation, which interrelates the proton-proton total cross section, elastic and inelastic scattering. 
The unitarity relation is, in fact, valid at each value of the impact parameter separately, that 
is 

2ImT;, = Y: KjT^,, (2) 

n 

where Tfi{s, h) is the transition amplitude to go from state i to state /. 

We follow a presentation by Pumplin after the original interpretation of Good and 



Walker ||15[. First we introduce states 0^ which diagonalize the diffractive part of the T matrix. 
Such eigenstates of diffraction only undergo elastic scattering. Let us denote the orthogonal 
matrix which diagonalizes Im T by C, so that 

ImT = CFC^ with (0fc|F|0j) = Fj 5jk- (3) 

Now consider the diffractive dissociation of an arbitrary incoming state 

\l) = Y C^k\<l>k). (4) 
k 

The elastic scattering amplitude for this state satisfies 

{t\lmT\t) = Y \C^k\'Fk = (F), (5) 

k 

where Fk = {4>k\F\(j)k) and where the brackets of (F) mean the average of F over the initial 
probability distribution of diffractive eigenstates. After the diffractive scattering described by 
Tfi, the final state |/) will, in general, be a different superposition of eigenstates than those of 



2 



\i) shown in (|). Suppose for simplicity, we neglect the real parts of the diffractive amplitudes, 
then 



dat, 



d(Je\ 



2lm{t\T\t) = 2 J2 \Cik\^Fk = 2(F) 

k 

mm\" = (e \Qk\'F,^ = {Ff (6) 



d% 
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It follows that the cross section for the single diffractive dissociation of a proton, 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 

Note that if all the components (pk of the incoming diffractive state \i) were absorbed 
equally then the diffracted superposition would be proportional to the incident one and again 
the inelastic diffraction would be zero. Thus if, at very high energies, the amplitudes Fk at 
small impact parameters are equal to the black disk limit, Fk = 1, then diffractive production 
will be equal to zero in this impact parameter domain and so will only occur in the peripheral 
b region. This behaviour has already occurred in pp (and pp) interactions at Tevatron energies. 
On the other hand, if there are, say, two diffractive channels with different eigenvalues, then 
the amount of inelastic diffraction increases with the spacing of the two eigenvalues. 



For instance, consider just two diffractive channels |jT8|, |T2|, ^ (say, p,N*), and assume, for 
simplicity, that the elastic scattering amplitudes for the two channels are equal. Then the T 
matrix has the form 

Im T = 1 - e"^/^ (8) 
where the eikonal matrix Q has elements 

nf,-, = VtoJ'ujf^i,. (9) 

The individual lo matrices, which correspond to transitions from the two incoming hadrons, 
each have the form 

The parameter 7(5, h) determines the ratio of the inelastic to elastic transitions. The overall 
coupling VLq is also a function of the energy and the impact parameter h. 

With the above form of cj, the diffractive eigenstates are 

= 7f + 1^*)) ' l'^^) = ^^\p)- ■ (11) 



In this basis, the eikonal has the diagonal form 




f2o r 



,mn 



(12) 



where m,n = (f)i,(j)2 and 



r = 



1 + 7 






1-7 



) 



(13) 



In the case where 7 is close to unity, 7 = 1 — e, one of the eigenvalues is small. 

3 Parton configurations of the difFractive eigenstates 

The simple two-channel model of Section 2 allows the prominent features of hard diffractive 
processes to be explained, which are beyond the scope of the single channel eikonal. First we 
note that the parameter 7, which determines the ratio of inelastic to elastic transitions, needs 
to be in the range 0.4-0.6 to be in accord with the experimental data on diffractive dissociation 
at moderate energies. Thus we know that there will be a big difference (1 ±7) in the absorptive 
cross sections for scattering in the two diffractive eigenstates. To be specific, in this work we 
use the results of the detailed analysis of the elastic and soft diffractive data that was presented 
in Refs. |rU|. There 7 was taken to be 0.4. 

In QCD the diagonal states correspond to quark and gluon configurations with different 
transverse coordinate^. For small transverse size r such (colourless) configurations interact as 
small colour dipoles with total interaction cross sections ~ r^. Thus, to a rough approximation, 
we can separate all the parton configurations of the colliding hadrons into those with small size 
and those with large size. In our two-channel example above these would correspond to the 
states 102) and respectively. 

It is informative to discuss the phenomenon in terms of the usual Reggeon diagrams. Assume 
that some "hard" diffractively produced stateQ "/i" is strongly coupled to state |02) and weakly 
to It follows from ([TT| ) that the Pomeron couplings of h to p and A^* satisfy 



and that the p and A^* intermediate states for double-Pomeron exchange contribution (Fig. 1) 
interfere destructively, since 



The cancellation which occurs for 7 ~ 1, happens, in this simple model, for all multiple-Pomeron 
exchanges. This phenomenon of "colour transparency" for small-size configurations has been 
known for a long time [pl|] . 

•^Partonic models of difFraction were originally introduced in Refs. pO| . 

^The state h should really be regarded as a third difFractive channel, but such a new state with a small 
production cross section gives a negligible contribution to n. 
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(14) 



(15) 
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Figure 1: The double- Pomeron exchange contribution to diffractive h production in the simple 
two-channel model of fP^. 



In order to specify the diffractive eigenstates and |02) we shall consider two simple 
models (A and B). It is natural to identify the component with the smaller absorptive cross 
section (that is |02)) with the state which contains less partons and which has a large typical 
momentum fraction x for each parton. From the QCD viewpoint, the small size component 
of the proton (where all the valence quarks are close together) has the smallest absorptive 
cross section, due to colour transparency. From the Regge viewpoint, the component with the 
largest absorptive cross section corresponds to the eigenstate (|0i)) with a larger number of 
partons in the small x region. Thus from both viewpoints we expect the component with the 
smaller cross section (smaller transverse size) to have a larger average x of each parton. At the 
moment, it is impossible to be more specific, and so to make numerical estimates we consider 
two alternatives. 

First, in model A, we identify the valence quarks with |02) with the smaller absorption, 
and the gluons and sea quarks with \4>i). Of course the model is oversimplified. It is clear 
that there is a part of the valence component with large size, while on the other hand the 
gluons and sea quarks contribute to the small size component. In general, one can write each 
partonic distribution /i(a;, Q^) {i = valence, sea, glue) as the sum of a small (S) and large (L) 
size component 

Mx,Q') = ff{x,Q') + //'(x,Q2). (16) 

In a model, where the probabilities of the S and L components in the proton are equal, as in 
Section 2, these components should satisfy the following sum rules, 

f'dxf^ix^Q') = I' dxf^{x,Q^) = 1 (17) 
JO Jo I 

[' dxxJ2 ff{x,Q') = f dxxY^fK^.Q^) = \. (18) 

JO ■ JO ■ / 

which follow from the conservation of valence quark number and energy respectively. 

We can therefore introduce an alternative model in terms of modified parton distributions 

/f'^(x, Q^) = P^'\x, Q') U{x, Q'), (19) 
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where the projection operators have the simple forms 

P-^ = (1 - xf^(Q'\ pf = I - i^^. (20) 

We determine the values of nj in order to satisfy the sum rules of ( p!?]) and (|18|). We call this 
model B. It turns out that both models A and B give rather similar predictions. We study the 
implications of the models in Section 5. 



4 DifFractive dijet production — a first look 



Recently CDF have measured diffractive dijet production for events with a leading antiproton 
at the Tevatron These observations, coupled with the diffractive measurements by HI 



22| and ZEUS at HERA, offer the opportunity to explore the diffractive framework in 
some detail. The processes are shown schematically in Fig. 2, in the absence of rescattering 
corrections. The lower parts of the diagrams, shown as Pomeron exchange, are to be understood 
as including multiple Pomeron contributions. 



TEVATRON HERA 




Figure 2: Schematic diagrams for diffractive dijet production at the Tevatron and for diffractive 
deep inelastic scattering at HERA. The rescattering corrections are omitted in these diagrams. 



If we ignore rescattering corrections, for the moment, then the cross section for diffractive 
dijet production of Fig. 2(a), integrated over t, may be written as 

^ = E/ FjpiOfriP,E^)f,ix,,E^)adpdx,d^, (21) 

i,k 

where a is the cross section to produce dijets from partons carrying longitudinal momentum 
fractions Xi and (3 of the proton and Pomeron respectively. This would correspond to the 
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Ingelman-Schlein conjecture |2||. Information on the diffractive structure functions ff{(3,Q'^) 



is obtained from measurements of the process of Fig. 2(b) at HERA F]p{^) is the flux 

factor taken for the Pomeron 



FAO = J dt ^5^, (22) 

where ^ is the fractional momentum loss of the recoil antiproton. In Regge theory, the coupling 
satisfies Cp = {g^iy/lGir, such that the total pp cross section is given by 



^tot(pp) = {gZr (s/so)^, (23) 

where sq = 1 GeV^ and A = ajp(O) — 1. When ^ is not too small the contribution of secondary 
Reggeons must be added. 

CDF present measurements of the ratio of dijet production for ii^T(jetl, jet2) > 7 GeV 
with a rapidity gap to that without a gap as a function of x = (3^ (the fractional longitudinal 
momentum of the p carried by the parton), for six ^ bins in the range 0.035 < ^ < 0.095 
with |t| < 1 GeV^ [|T7|. In the ratio, the terms fj{xi,E^)a cancel, assuming that single gluon 
i-channel exchange dominates the hard subprocess. Hence the data determine the diffractive 
structure function of the antiprotonQ 

^? = , 1 , r^'^ d^ FjpiO f3 [/f (/3, E|) + I /f (/?, E^^) 

(24) 

+ secondary Reggeon contributions. 

The CDF measurements of Fjj are shown by the data points in Fig. 3, together with five 
curves representing predictions of Fj^ based on various sets of diffractive structure functions, 
themselves obtained by fitting to HERA diffractive data. The structure functions ff {(3, Q"^) 
are evaluated at = 75 GeV^, which approximately corresponds to the average E^ of the 
CDF data. 

The prediction labelled by HI is obtained from the HI diffractive data, and corresponds 



to fit 2 of the HI collaboration [22]. The curve labelled by ZEUS (Pom) corresponds to the 
prediction obtained from ZEUS data in Ref. [p^p| - K does not include the contribution of 
secondary Reggeons. Note that the ZEUS data are in the region of very small ^ and thus are 
practically insensitive to these contributions. The curve labelled ZEUS' includes the secondary 
Reggeon contribution as determined by HI collaboration^. A comparison of the two latter 
curves shows that the non-Pomeron "background" is rather important in the ^ region covered 
by CDF (about 50% of the total contribution). These three predictions are representative of 



^Here we define the Pomeron flux slightly differently to Ref. by including Cjp in (22). 

^Note that the curve in Fig. 3 for the ZEUS structure function differs from that calculated in p5| |. 

'^This procedure may not be completely consistent as values of the Pomeron intercept are different in the 
analyses of the HI and ZEUS data (see Refs. This can lead to a modification of the secondary Reggeon 

contribution for the ZEUS parametrization. 
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those obtained from the various sets of diffractive structure functions that are available 126 



They illustrate the large uncertainties in the predictions of the shape of Fj^ at large P, and in the 
overall normalisation. On the other hand, the shape predicted for /3 < 0.15 is well determined 
to be with S = 0.4 — 0.5, and differs markedly from the measured 5 ~ 1 behaviour of the 
CDF data. 




0.04 



0.05 0.08 0.1 



P 



0.4 



0.6 



Figure 3: A comparison of the measured CDF dijet diffractive distribution as a function of 
f3, with different predictions obtained from analyses of HERA diffractive data assuming Regge 
factorization and that rescattering corrections are neglected. The shaded region on the CDF 
data shows the band of uncertainty shown in Ref. W7 . 



The diffractive gluon distribution is the main contributor to the predictions of Fjj. Although 
the diffractive quark distributions are well-measured at HERA (since the photon couples directly 
to the quark), the gluon distribution is determined from the detailed behaviour of the 
experimental data using QCD evolution. Moreover, the uncertainties in the diffractive structure 
functions ff{P, Q^) are amplified by sizeable differences between HI and ZEUS diffractive 
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data in certain /?, domains. These uncertainties mainly affect the gluon distribution and are 
responsible for the ambiguity in the predictions for the shape of F^- at large /3, and in the overall 
normalisation. All the predictions give similar shapes for < 0.15, because they are determined 
mainly by the QCD evolution. We will therefore study this evident difference between the CDF 
and HERA shape of at small as well as the difference in overall normalisation]^, which 
are clearly not reproduced in the naive model based on Fig. 2. 

Although we see from the curves shown in Fig. 3 that, at present, there are large uncer- 
tainties in the Pomeron structure function measured at HERA, curves I and II should provide 
a realistic illustration of the range of acceptable values, for the reasons given above. These 
two alternative curves are used in the analysis of the CDF data presented below. The curve 
I corresponds to the function 4.6(1 — and is very close to the parametrization of 



Capella et al. (with secondary reggeons taken from |22|), while the curve II corresponds 
to the parametrization 2.5(1 — which we choose to account for possible variations 

due to the uncertainty in the secondary Reggeon contribution. Recall both parametrizations 
corresponding to = 75 GeV^. 



The discrepancies between the Tevatron and HERA data were discussed in where it 
was emphasized that the survival probability of the gap is (i) small, and (ii) dependent on 
the value of /3. Physical arguments were presented which qualitatively reproduce the scale of 
normalisation and some trends in the /? dependence at large /3. Note that the effects causing the 



observed /5-dependence of the diffractive structure function considered in |2^ and in this paper 
concern different regions in /3 and are of different dynamical origin. While the fall-off at /? — 1 
is attributed in dominantly to Sudakov suppression effects, in this paper the variation of 



the shape of the /3-distribution is explained mainly by the competition between the different 
parton configurations. In this way, we present below a two-channel model prediction, based on 
0, which turns out to be in surprising agreement with both the normalisation difference and 
in the shape of the distributions at low /3. 



5 Diffractive dijet production including rescattering 
effects 

To explain the main features of the CDF diffractive dijet data it is sufficient to consider the two- 
component diffractive models introduced in Section 3. In model A we assume that the sea quarks 
and gluons mainly occur in large-size configurations of the incident proton, while the valence 
quarks occupy predominantly small-size configurations. This is, of course, an oversimplification 
of the real situation, but we find even this simple physical model is able to account for the 
behaviour of the data. 

^Note that earlier CDF results ||2^ on diffractive W boson, dijet, &-quark and J/V" production rates, using 
forward rapidity gap tagging, have already provided evidence against approaches which do not account for 
rescattering effects. 
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The two-channel generahsation of (|ID gives, using modelQ A of Section 3, the survival prob- 
ability of the gap^ 









2 g-nsca(s,b)^ 


I d%{\M.\ 


|2 + 


lA^sca 





where A^v,sea are the probability amplitudes (in impact parameter space) of the hard diffractive 
process corresponding to the valence quark and to the sea quarks and gluons respectively. The 
functions f2j can be parametrized in the form^ 

= ^'/'"^^ e-^V4B^ (26) 

Ann 

with i = V, sea, and where the slope of the Pomeron amphtude is 

B = ^Bo + a' Hs/so), (27) 

with So = 1 GeV^. We take = 1 — 7 and K^ea. = 1 + 7, consistent with the simple physical 
model introduced above. The values of the other parameters were determined in a two-channel 
global description of the total, differential elastic and soft diffraction cross sections 0, in which 
the parameter 7 was fixed to be 0.4. 

First we indicate why the soft rescattering effects {Qi 7^ 0) of the model based on ( pSj ) modify 
the (3 distribution of the dijet process in a characteristic way. Note that the CDF measurements 
cover a narrow ^ interval, 0.035 < ^ < 0.095, and hence that the invariant mass squared of the 
diffractively produced state, = ^s, remains close to the average value 2 x 10^ GeV^ Also 
the mass squared of the produced dijet system, 

Mj = xif3M^, (28) 

see Fig. 2, does not change much compared to its average value of about 1 x 10^ GeV^ calculated 
for the CDF kinematical range. Thus Xif3 ~ 0.005 and so for j3 > 0.25 we have Xi < 0.02, 
whereas for f3 ~ 0.025 we have xi ~ 0.2. Therefore for large /3 (small xi) sea quarks and gluons 
will give the dominant contribution, while for small (3 the valence quarks play an important 
role. Hence the survival probability should increase as Xi increases and f3 decreases. 



^In model B the subscripts 'v' and 'sea' correspond to the components with the smaUer and larger absorption 
cross sections respectively. 

"'^"In fact in the calculations a more accurate formula is used which takes into account the inelastic rescatterings 
of both of the colliding protons, see Appendix B of Ref. |^. 

"'^"'^We show formula (|2^) in order to again simplify the discussion. In this simplified form the values of the 



parameters would be about (g^)^ = 25 mb, Bq ~ 8 GeV^^, A = 0.1 and a' = 0.15 GeV^^. However, in 
practice we use the more realistic VLi(i = v,sea) that were determined in the global description of total, elastic 
and soft diffractive data in the ISR to Tevatron energy range . In addition the pion-loop contribution in the 
Pomeron was included (that is the nearest t-channel singularity) , as well as the contribution coming from large 
mass single- and double-diffractive dissociation. These refinements are not crucial for the effects that we discuss 
here. 
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Figure 4: The predictions for diffractive dijet production at the Tevatron, obtained from two 
alternative sets of 'HERA' diffractive parton distributions I and II (of Fig. 3), compared with 
the CDF data [0. The upper two curves correspond to the neglect of rescattering corrections, 
whereas the lower four curves show the effect of including these corrections using model A (con- 
tinuous curves) and model B (dashed curves) for the diffractive eigenstates (|0i) of Sections 2 
and 3). 



The calculation of the diffraction dijet rate, incorporating the rescattering effects of (p5D , 
confirms these expectations, as shown by the lower pair of continuous curves (I and II) in 
Fig. 4. These curves are parameter-free predictions of the diffractive dijet rate based on the 
two-channel eikonal model of Ref. ||^ and on the diffractive distributions obtained from HERA 
data. The two models (A and B of Section 3) for the diffractive eigenstates (|0i) and |02)) 
give similar predictions to each other, as shown respectively by the continuous and dashed 
curves in the lower part of Fig. 3. We see that the pair of curves II satisfactorily reproduce 
the normalisation and the experimentally observed shape of the (3 distribution. Curves I also 
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give a satisfactory description at low /?; the difference at larger (3 just reflects the uncertainty 
in the 'HERA' diffractive distributions. Recall that the predicted shapes show an anomalously 
strong increase[^ 1//?^ with 5 ~ 0.8 — 0.9, for small (3, as compared with the 5 ~ 0.4 — 0.5 
behaviour given by the partonic distributions of the Pomeron. A possible change of -ErOet), 
due to a variation ^Et of the transverse energy of the underlying event with /3, was taken into 
account in our calculations. We took = C{1 — PY with C = 0.76 GeV chosen so as to 
satisfy the observed (AEt) = 0.54 GeV The origin of such /3-behaviour can be traced 

to the fragmentation of the gluon jet. It leads to a small ~ 10% decrease of the theoretical 
predictions for /? > 0.2. 

The overall normalisation of the prediction for the CDF dijet data, which is reproduced 
by the average value of the survival probability (pSj), is sensitive to the impact parameter dis- 
tributions, Aii{s,b), of the hard diffractive process. Such a comparison can therefore provide 
information on these distributions which, in turn, reveal the spatial structure of the hard pro- 
cess. Our curves are obtained under the same assumptions for the single diffractive production 
of a massive hadronic state as were used in Ref. That is, as for the minimum bias single 
diffractive process, but without the term a'ln(M^/so), since a' ^ in (LO) DGLAP evolution 
to the scale fx^ ~ 75 GeV^ of the hard subprocess. 

Our calculation of diffractive dijet production illustrates a crucial ingredient necessary in 
the description of rapidity gap processes. Namely that the survival probability of a gap can 
depend on xi of the partons in the proton (see Fig. 2(a)). This leads to many experimental 
consequences for processes with rapidity gaps. For instance, if diffractive dijet production were 
measured at higher (LHC) energies with the same jet threshold (Ej,), then the values of Xi of 
the partons from the proton will be much smaller throughout the same interval of (3. Thus the 
variation of jS'p with (3 will disappear, and the shape of the (3 distribution in this interval will 
be close to that measured at HERA. The effect that is observed at the Tevatron is predicted 



to occur at the LHC, but at much smaller values of f3, see (p8 



5.1 Other (3 dependent effects 

We also estimated other possible mechanisms that may influence the predictions of the dijet 
(3 distribution shown in Fig. 4. We discuss these mechanisms in turn below. None of them 
is expected to be significant, and anyway influences mainly the region of f3 ~ 1. Most of the 
effects tend to make the (3 distributions steeper and to improve the agreement with experi- 
ment. However we have not included them in our predictions so as not to obscure the main 
phenomenon discussed in our paper. 

(a) The mechanism shown by the diagram of Fig. 5(a) describes the situation where the 
Pomeron couples to the 'ladder' in the upper part of the diagram, rather than to the 

^^This increase is still somewhat weaker than that seen in the data {6 w 1) |l7). However the CDF data 
include up to 4 jets, while the theoretical predictions are given for 2 jet production. If the data are restricted 
to two jet production then the increase is less steep (and given by the lower part of the shaded band) and, 
in fact, in agreement with our /3 dependence. 
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(a) (b) 

Figure 5: (a) Another possible contribution relevant to the calculation of the survival prob- 
ability of the rapidity gap in diffractive dijet production, where the 'rescattering' Pomeron 
couples to the upper 'ladder', rather than to the incoming proton, (b) The contribution to 
the diffractive dijet cross section coming from the interference of diagram (a) with the Born 
diagram of Fig. 2(a). 



proton as we have considered so far. It may influence the value of the survival probability 
at very small Xi. The interference of Fig. 5(a) with the Born diagram of Fig. 2(a) leads 
to a contribution to the cross section shown in Fig. 5(b). It contributes when the rapidity 
intervals yi (shown in Fig. 5(b)) are large; note that yi + y2 = In l/xi. This effect is small 
at the Tevatron because of the lack of phase space {yi + 2/2^5), but it should be taken 
into account at the LHC. 

(b) Another possible non-factorizable contribution is where a soft gluon from the Pomeron 
couples to the upper partons or spectator quarks of the proton 0, Their dominant 
contribution may be summed and absorbed in the reggeization of the gluon, that is BFKL 
effects in parton evolution. The remaining contribution is strongly suppressed since, when 
the soft t-channel gluon crosses an s-channel parton, it changes the colour structure of 
the corresponding splitting kernel. For the singlet (gluon ladder) Nc is replaced by Nc/2, 
whereas for the non-singlet (quark ladder) Cp = {N^ — l)/2Nc is replaced by —l/2Nc. 
Using the double log approximation we estimate this effect increases the prediction by 
less than 10% for small P, and less than 4% for (3 > 0.3. 
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(c) For large [3 there is an additional Sudakov-like suppression due to QCD radiation from 



high Et jets |2^. Again the dominant contribution is already included in the effective 
structure of the Pomeron measured in DIS at HERA. However conventional DGLAP 
evolution does not account for double logs of the type {as/^Ti) ln^(l — which sum up 
to exp(— (05/471) C ln^(l — /?)). The effect may be important when /5 — > 1, but is negligible 
in the /3 < 0.5 domain of present interest. 



(d) Hadronization may change the longitudinal momentum of the high Et jet |24[. The 
general effect is to shift the hadronic jet towards the centre-of-mass of the hadronic state, 

in Fig. 2(a), and to reduce the effective value of (3. A conservative estimate is that 
the prediction for F^^ is changed by less than 10% for [3 < 0.3, although it works in the 
desired direction to steepen the (3 dependence. 

(e) The predictions depend on the spatial size of the triple-Pomeron vertex. The correspond- 
ing slope hp is small, but not well known. We use the same slope hp = 1 GeV~^ as in 
1^, but 6p = or 2 GeV~^ are not excluded. Moreover we may expect hp to be smaller 
for larger /3 — > 1, when the Pomeron couples just to the hard sub-process. If we take 
hp = [1 — (3) GeV~^ then the prediction is unchanged in the small (3 region, although it 
decreases by about 10% at /5 = 1/3. The sensitivity to the radius of the triple-Pomeron 
vertex indicates the importance of the experimental study of diffractive dissociation pro- 
cesses to better determine hp. 



6 Predictions for other hard diffractive processes 

The observation that the suppression factors can depend on the values of the momentum 
fractions Xj, carried by the partons in the colliding hadrons, has implications for hard diffractive- 



like processes in general. For example, diffractive VT-production at the Tevatron |31[] is mediated 
dominantly by valence quarks in the proton, and hence the survival probability for such a process 
is comparatively large, S"^ ~ 0.2 — 0.3. On the other hand, for diffractive processes mediated 
by the gluonic components of the colliding hadrons (such as 66, J/ipjip' or T production) the 
survival probabilities should be smaller S"^ ~ 0.06 — 0.1. 

An interesting application is to the production of two high pt jets (p^j. ~ ~P2t) separated 
by a large rapidity gap, as measured by both the DO and CDF collaborations at the 
Tevatron at two energies, ^/s = 630 and 1800 GeV. Both the quark and gluon components of 
the proton contribute in this case. However, in our approach the suppression factor depends 
strongly on the type of parton (model A) or on the x value of the parton (model B). For a 
fixed energy ^/s, the ratio of the quark to the gluon component increases as Et of the jets 
increases, and as the rapidity interval Arj between the jets increases. The relative importance 
of the quark component also increases as the energy ^/s decreases, simply due to kinematics. 
These features of the simple two-channel model give effects which move in the right direction to 
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explain outstanding puzzles in the interpretation of the DO and CDF data for jets separated by 
a rapidity gap Q . In particular, they help to understand the ^/s, Et and A?7 dependences 
of the colour-singlet (rapidity gap) fraction measured at the Tevatron ^ . Our model gives 



a natural explanation of the observation by the DO collaboration, that the suppression factor 
depends mainly on the x of the partons and increases strongly with x, see Fig. 4(d) of W^ - 

All of the effects discussed above can be studied in hard diffractive processes in p-nucleus 
collisions at RHIC and LHC Investigation of the A-dependence can provide new information on 
the strength of shadowing effects. Indeed for weak shadowing, the cross sections for coherent 
diffraction dissociation of a proton on nuclei behave as ~ A^/^, while incoherent diffraction 
cross sections behave as ~ A. In the opposite limit of very strong shadowing both cross 
sections have much weaker dependence on A, of the form ~ A^^^. Thus there is a strong change 
in ^-dependence of diffractive production on nuclei depending on the strength of the shadowing 
effects. 

We also note that the survival probability for central Higgs production by WW fusion, with 
large rapidity gaps on either side, is enhanced in the two-channel model in comparison with 



previous estimates [y, |TT|, which also included allowance for the survival probability. Thus 
in Ref. for WW H process at the LHC S"^ was found to be 0.15, while the approach 
of this paper gives 5*^ = 0.24. This is an important process because it appears that large 
Higgs configurations can be chosen such as to identify the Higgs over the possible background 
processes at the LHC ||36|, |3^ . The same survival probability is applicable to central Z boson 
production with a rapidity gap on either side, originating from t-channel gauge boson exchange. 
Therefore Z production at the LHC can be used to directly measure the survival probability 
of rapidity gaps relevant to Higgs production by WW fusion H]. 



7 Conclusions 

For hard processes with large rapidity gaps, we have demonstrated that the survival probability 
of the gaps has a much richer structure than is given by the simple one-channel eikonal approx- 
imation of (|l]). We introduced two-channel eikonal models in which either the valence quark 
and the sea quark (+ gluon) components of the proton have substantially different total cross 
sections of absorption a\°^{s), which we called model A, or alternatively, model B, in which the 
small and large size diffractive components are specified according to sum rules (|T^) and (|TBp. 
The two models give similar results, and predict that the survival probability of the rapidity 
gap has a characteristic dependence on the kinematics of the process. Data for diffractive dijet 
production at the Tevatron [0 enabled this kinematic dependence to be checked. Taking the 
parameters of the two-channel models which were previously constrained in a global description 
of total, elastic and soft diffraction data, we calculated the (3 distribution of diffractive dijet 
production. The results are shown by the lower four curves in Fig. 4. We see that there is gen- 



eral agreement between the predictions and the CDF measurements both in normalisation 
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and shape. In fact, the use of Pomeron structure function II quantitatively reproduces the CDF 
data |T^. We emphasize that if the Pomeron structure function were known unambiguously 
then we would have an essentially unique prediction for the Tevatron data, demonstrated by 
the small difference between the predictions of models A and B. 

Unfortunately, the agreement between the CDF diffractive dijet data and our calculations 
can be taken as a strong support for the low value of the survival probability 5*^, which leads to 
the rather pessimistic expectations for the missing-mass Higgs search at the Tevatron p, |37 |. 

As precise data for other hard processes with rapidity gaps become available, it will be 
possible to refine the model and to identify the parton content of the diffractive eigenchannels. 
We already showed that the simple two-channel model gave rescattering corrections which 
moved in the right direction to resolve discrepancies between the predictions and the data for 
processes which have so far been measured. In this way, as precise data become available, it 
will be possible to perform a quantitative study to (i) determine the partonic content of the 
Pomeron, (ii) check the QCD evolution of the Pomeron structure functions, (iii) confirm the 
universality of the partonic decomposition, (iv) determine cr*°*(s) for the different diffractive 
eigenchannels, and (v) measure the impact parameter distributions of the 'Born' amplitudes of 
the hard processes with rapidity gaps. 
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